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Abstract 

■ We resum the recently calculated second order kernel of the BFKL 

i equation. That kernel can be viewed as the sum of a conformally 

I invariant part and a running coupling part. The conformally invariant 

' part leads to a corrected BFKL intercept as found earlier. The running 

, coupling part of the kernel leads to a non-Regge term in the energy 

^ I dependence of high energy hard scattering, as well as a Q^— dependent 

0> ' intercept. 



1 Introduction 



Recently, the calculation of the second order kernel for BFKL|]1], |^ evolution 
has been completed in two (partially) independent calculations [0, The 
^ ■ second order kernel strongly modifies the leading order result for the BFKL 

intercept at all but the smallest QCD couplings. While this large intercept 
correction does not significantly alter our picture of BFKL dynamics, it sug- 
gests that even higher order corrections are important leaving us without a 
reliable theoretical calculation of the intercept which can, however, in prin- 
ciple be determined phenomenologically. 

Of course, the corrected BFKL intercept governs single-scale high energy 
hard scattering only after a resummation has been done to bring that inter- 
cept into an exponential form. It is the purpose of the present paper to do a 
more complete resummation of the complete second order kernel for use at 
high energies. 

*This work is supported in part by the Department of Energy under GRANT DE- 
FG02-94ER40819. 
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In Sec. 2, we develop a formalism which allows one to iterate the complete 
second order kernel an arbitrary number of times. It turns out to be conve- 
nient to separate the second order kernel into a conformally invariant part 
and a running coupling part. The resummation of the conformally invariant 
part leads to the result given in Refs.3 and 4. In Sec. 3, after some calculation, 
we resum the running coupling part of the kernel, and it also exponentiates. 
The total result is given in (29) and (30). The running coupling part of the 
kernel has two effects. It determines the scale of the leading order intercept to 

4iVefa2 Q,^^gg/-j -^j^gj^g q g^]^^ qi Qj-Q lY^Q virtualities of the ghions which hook 

into the (compact) particles for which the scattering is being calculated, as 
given in (1). The running coupling part of the kernel also introduces a non- 
Regge behavior in the scattering, the ^[a{ap — 1)6]^F^ term in (30). Indeed, 
this non-Regge term is necessary in order to be able to view evolution in a 
rapidity interval Y as built out of separate evolutions over regions y and 
y — y as indicated in (32). Although this non-Regge term formally looks like 
a next-to-next-to-leading term, in the exponent, it does not appear possible 
to generate such a term from a third order kernel and so we believe the result 
given in (30) is solid. If one writes this term as ca^Y^, c 5 for 3 flavors 
so that this term is reasonably small so long as a is reasonably small and Y 
is not inordinately large. Parametrically, the Regge-typc behavior in BFKL 
evolution begins to break down when Y ~ a~^^^ a value where perturbation 
theory is still valid, where the next-to-leading conformal corrections to the 
intercept are (parametrically) small, but where unitarity corrections are ex- 
pected to already be large. We note that E. Levin has previously arrived at 
Y < Q.^^l'^ as a criterion for the validity of the usual Regge-BFKL picture of 
high energy scattering. (See Eq.82 of Ref.5.) 

In Sec. 4, we calculate the next-to-leading corrections to the anomalous di- 
mension and coefficient functions at large orders of perturbation theory. The 
largest corrections to the coefficient function are determined completely by 
the shift of the branch point in the angular momentum plane from uo = ap — 1 
to the position of the corrected BFKL intercept, while the largest corrections 
to the anomalous dimension function are determined by the appearance of a 
pole, = as indicated in (41). 

Finally, in Sec. 5, we investigate the limitations on using perturbation 
theory to calculate single-scale high energy hard scattering. We do this 
by examining the factorials which appear in perturbation theory and which 
determine the region where the QCD asymptotic expansion can be used. 
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Eq.46 confirms the result found in Ref.6 using a fixed coupling approach to 
the diffusion present in BFKL evolution. 

2 Using the next-to-leading kernel to all or- 
ders 

The total cross section for the scattering of two compact colorless particles 
(heavy onia) A and B can be written as 

^-(^) -/$/l^*^(3)*B(g')/;_rS(^r«»(3.9') m 

where we use the notation of Ref.3. G,^ obeys the equation 

uG^iq, q') = 5(q - q') + / d'qKiq, q)G,(q, q') (2) 

with 

X(q,qO=i^W(q,q')+i^(2)(q^q/) (3) 

where the superscripts in (3) indicate the leading and next-to-leading BFKL 

kernels, respectively. If G^-^-* represents the leading BFKL result, the result 
of solving (2) taking only K^^^ as the kernel of that equation, then G^j can 
be written as 

00 

G.(q,q') = EG^^-(q,q') (4) 

N=0 

with 

G^. = G«[ir(2)G«]^ (5) 

where all of the products in (5) are understood to be convolutions as indicated 
in (2). It is convenient to write G'^^ as 
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with x(A) = 2ip{l) — '?/'(A) — — A), and where we neglect the dependence 
on coscj) = (q ■ q'/qq') a dependence which disappears when a{s) is evaluated 
at large s. The A-integration in (6) runs parallel to the imaginary axis with 
ReX = 1/2. 

It is useful to express the kernel K^'^^ as a conformally invariant part and 
a "running coupling" part using the result 13] 



with 



and 



d'q'K^'\q, q')(^)''-' = k^'\q) = k^nf{\) + Kc{\ q) (7) 



KonM) = -i(^^^)2c(A)x(A) (8a) 

4 TT 



krciX, q) = -^^^^^l^bx{X)in q/fi {8b) 

TT 

where c(A) is defined in Ref.3, and where b = Define 

G^{q,Y,q')= J ^e-yCM^.q'). (9) 

along with similar definitions for G and G^^^ Then, substituting (6) in (5), 
using (9) as well as the identity 

duj e^^ 1 



duj -iX- duj, d^^^y^-y^-^) (.^{y-vn) 

with = yo<yi<y2^ ■ ■ ■ <2/Ar<^, one finds 

GN{q, Y, q') = J G^'\q, Y - y^, qN)d\NdyNk^^\qN) 

N-l 

n G^'\q,+,,y,+,-y,,q,)d%dy,k^'\qi)G^'\q^,y,,q') (11) 

i=l 

where, in the saddle-point approximation. 
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with D = ^^^iV,C(3) and ap - 1 = ^"^^J^^° £n2, where (12) is vahd for large 
y. Assuming 0<<yi<<y2<< • • • «yN«Y in (11) one gets 

Jap-1)Y I 

GN{q,Y,q') 



IW^VF^. ^""'1 Y - y/''^^ ADiy^-y^J'""^^ AD{Y - y^Y 

(13) 

where itj = In qij /i, Uq — £n q'/ /i and u — in q/ /i. Now 



1 1 

k^'^\ui) = Konfi^) + krc{^,Ui) (14) 



with 



k^4) = ^'^' c(i) (15a) 

Kc{^, Ui) = -2a{ap - 1)6 Ui. (156) 

It is convenient to sum all orders of kconf for a given order of krc- Thus, 
one may write 

«(.y.,o (16) 

27iqqWATi DY 

with 



'AT 



N 

n 



[47rD]^/2 y ^y, - y,_. 



' ^ exp\- ~ ] 



■exp[- 



■4L>(y-yjv)^ 



(17) 
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Choosing Y = C-n-^ and using (16) and (17) in (1) gives the result of summing 
the leading and next-to- leading kernels to all orders for the cross section aAB- 
The first factor on the right-hand side of (16) is the answer given in Refs.3 
and 4. 



3 Evaluating the running coupling contribu- 
tion 

In order to evaluate In one can introduce a factor exp{Yjf^i JiUi} in the 
integrand of (17), expressing each Ui coming from krc{ui) as ^ acting on the 
integral, and then setting the Jj = at the end of the calculation 



In = [-2a{ap - 1)6] 



N 



N 

n 



d 



dy^ , Y 
duii 



It is useful to make the change of variables 



(18) 



Mo 1 



^0 



yi^x{ui - uq) - yi{ui+x - uq) 1 



yiyi+i{yi+i - yi) v4L' 

u 1 



,z = l,2,---,iV (19) 



\/y VW 

where yN+i = YjUn+i = u. The Ui are expressed in terms of the z[s by 



Ui-uo = yi 



N 



I Vj+i-Vj ^ z-Zq 

yjVj+i Vy 



'AD 



u-uo^Y^^V4D. 



Vy 



(20a) 
(206) 



After rescahng J, \MD J, we have 
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IN 



N 



[-2a{ap - l)bfexp{ 



TT 



Il-gj, J dyidziexp{- 



Vi+l — Hi y-v J.^.j 
yi+lUi j=i 



.exp{^=£n^{Ji + • • • + Jjv) + -yi=^in^{Jm + ■■■ + ^Jvyiv)}L=o (21) 



'AD H V4:DY q' 

where, as always, the y-integration is ordered 0<yi<y2<---<yN^y and where 
it is understood that all the ^ factors are to be put to the left of all the 
Jj terms in the various exponents in (21). After performing the integration 
over dzi and after some simple algebra formula (21) becomes 



IN 



[-2a{ap - l)bfexp{ 



ADY 



^n2^}(4D)^/2 

q 



LLgjJ dytexpi- }^ JiJj } 

i=l * j=i+l 



.exp{-^£n-{Ji + • • • + Jjv) + -y^Jn^{Jm + ■■■ + JivZ/iv)}L=o- (22) 



'AD A* 



'ADY q' 



Moving the terms linear in J in the exponent to the left of the factors 
gives 



1 n ^ 

Ir, = [-2a{ap - l)hf exp{-^£n'^]{ADri' U 



1=1 



j dyi 



in'- + ^y^in^ + ^ 



AD A* VADY q' dJi 



exp{\ E 



Now consider the expression 



(23) 



Pk = iil dyigjexpi- 2^ JiJj 



i=l 



(24) 



Pk is non-zero only for even k in which case 



7 



p 1 1 /• , , E n ViiY - yj) 



(25) 



2t (f)! - ^^.-^.<, 

where the sum goes over all permutations of the pairs i and j. Eq.25 is readily 
evaluated as 



1 1 



2k/2 (|^ 



Y 



dy2 / dy 
Jo 



yi{Y- y2) 
Y 



- k/2 ^ 




k/2 


\ " (1)! 


2-4! 





(26) 



Rewriting (23) as 



7. = [-M<^P-m''exp{-^in^^}m'^'EPkj^;^ 

(27) 

we obtain 



q; 



(ap - l)6)2DF3|fc/2 



where the sum goes over even values of k. Using (28) in (16) gives 

(,E{q,Y,q') 



G{q, Y, q') 



27r qq'VInDY 



(28) 



(29) 



with 



Eiq,Y,q') = {ap{qq') - 1)(1 - ^c(i))F + |(a(ap - 1)6)^^^ - 

(30) 

where we have used 

{apiix^) - 1)(1 - a{ix^)h{uo + u)) = q;p(??') - 1. (31) 
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The —{ap — l)^;;;^c(|) term in (31) gives the next-to-leading correction 
to the BFKL intercept ||, |. This term comes from the conformally invariant 
part of the next-to-leading kernel, as given by (8a) and (15a). The running 
coupling part of the next-to-leading kernel gives the scale of the leading 
BFKL intercept, the ap{qq') — 1 = ^^in2a{qq') term in (30), and also the 
non-Regge term, the a^Y^ term in (30). The appearance of a term in the 
exponent signals a breakdown]^] of the usual picture of high energy scattering 
where powers of Y and the factorial denominators associated with them come 
only from a strongly ordered region of longitudinal phase space. The a^Y^ 
term in (30) is, however, purely perturbative and the overall coefficient is 
reasonably small, y[a(ap — l)6]^y^ 5a^Y^, so that this term is likely not 
too important for present phenomenology. We note that the expression (29), 
with E given in (30), obeys the consistency condition 

G(g, Y, q') = J d'kGiq, y, k)G{k, Y -y, q'), (32) 

and that the a^Y^ term in E is crucial for (32) to hold. Indeed, one can view 
(32) as requiring the a^Y^ term, exactly as given in (30), once one has arrived 
at the scale qq' as the appropriate scale for the leading BFKL intercept. 

4 The anomalous dimension and coefficient 
functions at large orders [7- 10] 

Certain higher order corrections to the gluon anomalous dimension have al- 
ready been calculated]^. Here we do a much simpler calculation. We shall 
calculate the dominant large order terms at next-to-leading level for both the 
coefficient function and the anomalous dimension function. These large order 
terms are insensitive to the scheme in which the energy scale is introduced. 

To determine the coefficient function C^{a) we take q = q' = fi in (29) 
and write 

nqq'G{q,Y,qX^^,^^ = /^C^e^^- (33) 
Keeping only the first next-to-leading order term one gets 
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2V47rDF 

where 



1 - '>'\ d)y) = / ^(C.„ + (34) 



^^0,. = E (^)^-C^r, ^^1. = E (^)^Cr^ (35) 
The Co^"* are given by 

C^r.~oo^(^)--^^i_ (36a) 
VN aNjir 7567rC(3) 



while a simple calculation, using (34), leads to 

C^N-^ - ^Cr''. (366) 

In order to determine the anomalous dimension function at next-to-leading 
level one can take q = fi,q' = Q and evaluate (29) to first order in in Q^///^. 
One finds 



_p(ap~l)Y I. J 1 

^^=_a(ap - l)hY^nQ/^, = J £-[^(7. - ^)]Jn QV/^'e^^ (37) 

where [Ca;(7w — |)]i indicates the next-to-leading corrections to the product 
Cuiilu,- I). That is 

[C.il.-hl = t[Cil-hr\^f. (38) 
Using (37) it is straightforward to determine 

which, using (35) and 
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along with||] 7,3 



70,a; = J2 7o" 
N=0 

(TV) _ 1 ^(TV) 



aN, 



c\TV 



, > 7i,< 
leads to 



1:7! 



(TV) 



aNr,aN, 



N=l 



TC TTLU 



(40) 



bn ap 



\N-1 



(41) 

The correction for c[^^ given in (36b) simply corresponds to a shift in 
the branch point in from = [16D{uj — {ap — 1))]~2 
to = {16D[cj-(ap-l)(l-^^c(i))]}-i The A^^/a factor in (41) indicates 
the start of the appearance of a non-Regge term. The fact that Ci^^ /Cq^^ 
and Ti^Vto^'' large |ri[| at large N does not in itself signal a breakdown 
of perturbation theory since these terms are known exactly and, at least in 
the case of 7!^'', lead to corrections which are not particularly large. 



5 Limitations on the use of perturbation 
theory 

Because of the diffusion inherent in BFKL dynamics it is clear that the 
perturbative approach to high energy single-scale short distance behavior 
must breakdown at sufficiently high energy. Estimates of the energies at 
which that breakdown occurs were given in Ref.6. We are now in a position 
to see this breakdown in a more detailed way than has been done 
previously. The idea is simple. Perturbation theory should itself indicate 
when it is breaking down by generating factorial terms which indicate 
where the asymptotic expansion, the perturbative expansion, is reliable. 
Suppose we take the N = 1 term in (5), but now instead of taking the 
running coupling correction to be given by (8b) we work to arbitrary order 
in running coupling corrections 

~krc{X,q) = ^^^x(A) f:(-26«(^2)£n g^)'^. (42) 

^ n=l 
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For simplicity we take q = q' = /j, in evaluating G. Thus, we replace Ii in 
(16), (18), and (21) with 

DC 

h = E (43) 
n=l 

where 



[-2Q;(/i2)6]"(ap- 1) 
V47rL> 



rY PCX 
Jo J-c 



dzie 



yiY 



(44) 



A simple calculation leads to 



h,n 



lap 



2 n + 1 



^^^''(4aVL>y)"/2r(-) 



(45) 



for n even, and Ji „ = for n odd. Of coTirsc the sum indicated in (43) does 
not converge, and this is not surprising. So long as 4:b^a^DY«l, Ii^^ is 
small for small values of n and the divergence of the sum in (43) is a 
standard renormalon problem. One can keep those terms in n so long as 
l^fi±!- < 1. with higher orders in n being discarded. However, when 

Aa^b^DY ~ 1 wc are not allowed to safely use any of the terms in the 
perturbation series given in (45), and thus the whole perturbative approach 
breaks down. The criterion for being able to use perturbation theory 
without too much contamination from infrared regions then is 



WDb^Y<l 



(46a) 



or 



Y<- 



TT 



14Ar^C(3)&'"Q;3(/x2) W 

exactly as found in Ref.6. // characterizes the scale of the single-scale hard 
process in question. 
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